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Staggered Domain Wall Fermions (SDWF) combine the attractive chiral properties of staggered fermions with 
those of domain wall fermions. SDWF describe four flavors with exact U(l)xU(l) flavor chiral symmetry. An 
extra lattice dimension is introduced and the full SU(4)xSU(4) flavor chiral symmetry is recovered as its size 
is increased. Here, the free theory of SDWF is described and a preliminary discussion of the interacting case is 
presented. SDWF may be well suited for numerical simulation of lattice QCD thermodynamics. 



1. Introduction 

Numerical simulations of QCD thermodynam- 
ics, using current methods and supercomputers 
do not realize the flavor chiral symmetries as well 
as one would like. More specifically the study of 
the relevant critical phenomena close to the fi- 
nite temperature transition may be affected. It is 
reasonable to investigate alternatives to current 
methods especially if they combine some of their 
attractive features. 

Staggered fermions, and their variations, de- 
scribe four flavors (in four dimensions) and have 
an exact U(l)xU(l) flavor chiral symmetry. The 
full SU(4)xSU(4) symmetry is recovered when 
the lattice spacing a approaches zero. Two fla- 
vor simulations of QCD are done by taking the 
"square root" of the staggered determinant. Be- 
cause the U(l)xU(l) chiral symmetry is exact at 
any a staggered fermions have been attractive for 
studies of spontaneous chiral symmetry breaking 
phenomena. On the other hand only one pion in 
the multiplet acts as a Goldstone boson due to 
the flavor breaking. Even the highly improved 
asqtad action may require > 8Tc to get all 
pions lighter than the kaon mass ||l| . 

A few years ago new methods were developed 
with good chiral symmetry properties: please see 
the review articles in [^-^ and references therein. 
One of them, domain wall fermions (DWF), uses 
standard Wilson fermions but introduces an extra 
dimension with Lg sites and free boundary con- 



ditions. As a result massless surface modes de- 
velop, with one chirality bound exponentially on 
one wall and the other on the other wall. At fi- 
nite Ls chiral symmetry is broken by their overlap 
but it is restored exponentially fast as Ls is in- 
creased. At Ls — oo (the overlap formalism |^,^ ) 
the chiral symmetry is exact, even at non zero a 
and topological zero mode effects are reproduced. 
These remarkable properties make DWF a good 
candidate for QCD thermodynamic studies. For 
such studies see 

EMI 



Both staggered and DWF are very close to re- 
alistic simulations of QCD thermodynamics. The 
present work is an attempt to combine their at- 
tractive features. SDWF are staggered fermions 
defined in space time with one extra direction and 
free boundary conditions. As with DWF, surface 
modes develop with the plus chirality of all flavors 
localized on one boundary (wall) and the other 
chirality on the other boundary. By properly 
identifying the various flavors on the boundaries 
one can construct four flavors of Dirac spinors 
with exact U(l)xU(l) at non-zero a and finite 
Ls- The full SU(4)xSU(4) symmetry is recovered 
in the Ls oo limit even for non zero a. 

Here, the free theory of SDWF is described and 
a preliminary discussion of the interacting case 
is presented. At this moment it is not clear to 
what degree the localization properties persist in 
the strong interaction regime where lattice QCD 
thermodynamic simulations are presently done. 



2. Staggered Domain Wall Fermions 

We start our construction of the free SDWF 
action by writing the free staggered action in mo- 
mentum space and in the hypercubic notation of 
Kluberg-Stern et al. Q 



Ssi = EfcQ(fc) W ® 1)+ Hk)] Q{k) 



(1) 



may be simultaneously diagonalized. This con- 
straint alone effectively restricts the allowed pro- 
jectors to the ones we have chosen. We define an 
additional set of operators to project out specific 
flavor components as well: P+± = {P+ (1 ± 
^5))/2. The solution is separable and 0(s) is the 
s-dependent part. We write it block notation us- 
ing the new projectors: 



^(fc) = EJisinA;^(7^®l) + 6^(75 ®C5m)] (2) (/-(sf = (^^+, ^/-T-' 



(6) 



with — {1 — cosk^)/2. The lattice coordi- 
nates Xf_i are decomposed into the coordinates of 
2'^ hypercubes and coordinates of sites within 
a given hypercube = (a:^ mod 2) S {0, 1} so 
that Xf_i = 2y^^ + A^^ . The momenta are recip- 
rocal to the hypercubes y^. 

To complete the construction of the free SDWF 
action, we introduce an extra dimension in the s 
direction with free boundaries, sum Ssi over the 
s-slices and add terms to the action that project 
on (75 (g) 1) chiralities 

W = Es [Ssi + EMk, s)D^Q{k, s)] 



^5(5, 



s+l,s' 



-l,s' 



(3) 
(4) 



Explicitly, the chiral projectors are P± = ((1 ± 
75) ® l)/2. Even though breaks the 

U(l)xU(l) symmetry generated by (75 (8)^5), we 
will show in section ^ that an extended symme- 
try exists to protect mass terms like m(l ® t-)5s,s' 
from additive renormalizations. Thus, we find it 
natural to set m — > — and drop these terms from 

as ^ 

the action. However, it is important to keep in 
mind that at higher energy scales, fermions prop- 
agate like massive staggered fermions in d dimen- 
sions. 

Fully chiral fiavor symmetric modes should ex- 
ist at lower energy scales in the SDWF theory if 
we can find normalizable states that satisfy the 
following zero mode equation 

{i;[P+^s+i,s' + P-5s-iA 
+ T.^,{l^®i^n^)b^.5s,s]^{k,s') = Q. (5) 

From this equation, it is easy to see that the P± 
projectors in the s-dependent part should com- 
mute with the flavor breaking part so that each 



where 4> h(*) = P v4'{s), etc. In this notation, 

we can write 



( t 

-fit 



V 



B 



St 



•(7) 



The solutions to (after iterating one time) are 



b±+{s±2) 
6±_(s±2) 



(8) 



Solving the equations relating nearest neighbor s 
sites is more complicated because the flavor com- 
ponents mix and will be discussed elsewhere . 

For free fermions, [-8,5'''] = and a^BB^ , 
a^B^B are both proportional to the identity with 
eigenvalue 

X{alBB^) ^ X{alB^B) 



(9) 



If we require that aiE/i^^ ^ ti^en for a semi- 
infinite s direction, s>0, only (t>+± is normaliz- 
able, while (f>-± is not. However, this is not 
enough to ensure that the doubler modes are 
not present. One must further require that the 
above condition excludes momenta with compo- 
nents larger or equal to tt. This is satisfied pro- 
vided that a5>l. If one wanted to further restrict 
the momenta to tt/2 then the requirement is that 
a5>2. 

Most of the recent work in the overlap formal- 
ism is related to taking the as— >0 limit. Clearly, 
additional terms must be added to our SDWF ac- 
tion to maintain the cutoff of doubler momenta as 
05— >0. Possible terms are currently under study. 
In any construction of a staggered overlap Hamil- 
tonian, it will be important to demonstrate that 
naive fermions are not recovered in the chiral limit 



and that the U(l)xU(l) symmetry remains exact. 
Of course, it is always possible to use the overlap 
formalism directly with finite 05. 

When adding interactions to a free staggered 
action, it is usually a good idea to first tran- 
scribe the free action in terms of single com- 
ponent per site fields in position space. The 
unitarily equivalent hypercubic basis of Daniel 
and Sheard is ideally suited for this pur- 
pose. In this basis, fermion bilinears are writ- 
ten XAiy)ilS £.F)ABXB{y), where the fermion 
fields XA{k) are directly the single component 
fields on the corners of the hypercube y: x(^) = 
xC^y + ^) = XA{y)- The components of the spin- 
flavor matrices are given by 

hs®^F)^s = ^tr[j\jslBlU G {0,±1} (10) 

where S and F are d-dimensional binary vectors 
like A and B and 75 = 7f ' x ■ • • x 7;^^* . 
The staggered P is the usual one 

^-E^(-i)'"^*"'[^x+a,.'-'5.,.'+a] (11) 

where ri^{x) — "^^^^Xi^. For the projection 
terms in D5 proportional to (75 (X" 1), we have 

X(y,s)(75 «) l)x(2/,s± 1) 

= {-ir^^^x{x,s)x{x,s±l) (12) 

where (p{x) = d/2 + X^^l^i X2ij,-i and x is the 
opposite corner of the hypercube: = + 1 — 
2{Xf, mod 2) V /i. 

3. Symmetries 

When constructing the SDWF action, it is im- 
portant to preserve all the symmetries of the 
massless staggered action |Q . Of course, adding 
any new terms to the staggered action will likely 
break some of those symmetries, so we have to 
find new symmetries that involve the extra di- 
mension. However, we will still have to show that 
these new symmetries in the space with an extra 
direction are maintained by the action that has 
all fields integrated out except the "light" ones at 
the boundaries. We present here the symmetry 
transformations in the hypercubic notation for 
the space with an extra dimension. 

Rotations by tt/2. These rotations are in planes 
perpendicular to the extra dimension and the 



transformations are the same as the original stag- 
gered ones. 

Ii-parity. These transformations reflect the d — 
1 spatial axes perpendicular to the spatial axis in 
the jji direction. D5 is not invariant under this 
symmetry unless we also reflect the s direction as 
well. If we define the reflection operator TZg^s' = 
Sl^^i-s.s' , the transformation is 

Qiy,s) (jt, ^5)'R's, s' Qiy,s'), 

Q{y,s) ^ Qiy,s') 7^,^,(7^ 0^5). ^ ' 

Shift by one lattice spacing. By inspection, one 
can see that (^ij) is invariant under shifts by one 
lattice spacing in any spatial direction /z. Be- 
cause of the additional structure imposed on the 
lattice by the hypercubic formulation, the sym- 
metry transformation, while still valid, is compli- 
cated. For SDWF, there is an added complication 
that some parts of the transformation require a 
reflection in the s direction (with s indices sup- 
pressed): 

Q{y) -> i [(1 ® e^) - (7,,5 ® ^5)Tl] Q{y) 
+ \ [(1 ® Cp) + (7^.5 ® Qiy + fi), 

Q{y) -> Q{y)h [(i ® - ® Cs)] 

+Q{y + [(1 ® 00 + 7^(75m ® ^5)] • (14) 

U(l)e xU(l)o chiral rotations. The residual chi- 
ral symmetry of staggered fermions involves mak- 
ing separate chiral rotations on even and odd 
sites. Terms in are not invariant under these 
rotations unless we extend the notion of even and 
odd, including the extra dimension. First, we de- 
flne the operator Ss,s' = {—^Y5s,s' and then the 
extended even/odd projection operators 

P, = |[(l®l)+5(75 0^5)], 

Po = i[(l®l)-5(75®^5)]- (15) 

Using these projection operators the chiral trans- 
formation is 

Q{y) ^ (e'»=Fe + e*^°Po) Q{y), 

Q{y) ^ Q{y) {e-''"Pe + e^'-Po) ■ (16) 

4. Flavors of SDWF 

From section || one sees that for a finite ex- 
tra direction with Lg sites the P_|_ components of 



all flavors are localized around s—0 while the P_ 
components are localized around s=Ls — l. From 
section however, we mentioned the importance 
of constructing an effective d-dimensional field q 
from the surface states consistent with all the 
SDWF symmetries, particularly the U(l)e xU(l)o 
chiral symmetry. For example, to project flavor 
components with P++, one should choose s near 
zero. If s=0 is chosen, P++q{y) = P++Q{y,0), 
then these components also belong to the Pe part 
of the fermion field. Therefore, to project fla- 
vor components with P ^ one is not only re- 
stricted to choose s near Lg — l but also choose s 
so these components belong to the Pq part of the 
fermion field. Then, components P±+q will not 
mix even for finite because of the even/odd 
symmetry. In this example, one would like to 
pick P^+q{y) = P_4.(5(j/, s) with s being even 
and near Lg — l. So, if Ls is odd then s=Ls — \ is 
a good choice. However, if Ls is even, then one 
should choose instead. 
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(17) 



Q 



PeQ 



Using the block notation of 
for even is sketched in (17f). 



PoQ 

@, an example 
In this equation 



Q{s—0) is at the top and Q{s=Ls — l) is at the 
bottom. The capital letters denote one of the 
"correct" choices. On the other hand if Lg is odd, 
say 3 then 
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(18) 



Q 



PeQ 



PoQ 



We note that other choices for selecting flavor 
components near the boundaries are certainly 
possible. We also note that the even Lg choice 
in ( p^ does not transform simply under the ft 
reflections of section ^ and requires modification 
to make it consistent with the other staggered 
symmetries. This will be discussed in more detail 
elsewhere |T^ . 

5. The SDWF propagator 

The SDWF propagator for the Dirac matrix 
given in section |^ is presented here in a general 
form. The detailed form will be presented else- 
where [|l^. A degenerate mass term that explic- 
itly mixes chiralities is added. For odd Ls it has 
the form 



Ms,s' = -imfSsx,-iSo,s'P+ 
+imfdsfi5L,-i^s'P-- 



(19) 



For even Ls it has a form according to the discus- 
sion in section ^. The propagator in momentum 
space has the general form: 

D-\s,s')= (20) 
[Gi + mfG2]e{s - s') + G3e{s - s' - 1) 



where Gi, G2 and G3 are functions of momenta 
and are proportional to the identity in their flavor 
indices. Also, Gi anti-commutes with (75 ® U) 
while G2 and G3 commute with (75 (g) 1). The 
flavor mixing is in the function e: 



e{x) = (1 (8) 1), {x even) 



{x odd) 



where |6| — \/J2fj. ^Ji- For s — s' even and 



(21) 



=0 



the propagator anti-commutes with (75 (>D 1) and 
has no flavor mixing except for the last term in 
(pO|). In this term e(odd) breaks flavor in ex- 
actly the same way as free staggered fermions. An 
exact U(l)xU(f) symmetry is maintained. The 
matrix coefficient G3 vanishes exponentially fast 
with Ls for s, s' near opposing boundaries and 
therefore as — * 00 with m/=0 the propagator 
anti-commutes with (75 (g) 1) and has no flavor 
mixing provided s — s' is even. This is in accor- 
dance with the discussion in section ^. 

Finally, if we add to 6^ a constant term l/as — 
TTT'O- bi_i = {1 — cosk^ -|- l/as — "io)/2, the effective 
mass meff has exactly the same form as in Wilson 
DWF but the decay coefficient is now in terms of 
|6| instead of 6 = ^^[^ — cos{pfj,)] -f 1/05 — mo. In 
this case, the localization condition for mg is the 
same as in DWF and the a^^Q limit can be taken 
the same way as in DWF. Furthermore, this term 
may be needed to cancel any renormalization of 
the flavor breaking term. This term would appear 
in the action as (l/os -mo)/2 Q{l5(E)^5tJ.)Q, a 
term that is not invariant under shifts by a single 
lattice spacing. However, such a term, or a gen- 
eralization of it, might be needed for the above 
reasons. 



6. The SDWF transfer matrix 

We can use the technique of Neuberger Q to 
rewrite the free SDWF determinant in a form that 
allows us to quickly identify the transfer matrix. 
We use the same basis as in (|7|) . After interchang- 
ing various rows and columns of the SDWF ma- 
trix, the determinant is equivalent to the deter- 



minant of the matrix 



Pi 



V 



(22) 



where all of the and Ps are the block triangular 
matrices 



as 



B 

-C 





l/as 



, Ps 



l/as 




-B 



.(23) 



For aL^-i and Pq, l/as is replaced with — /x/as 
so /i is a parameter that controls the boundary 
conditions: /i = ±1 for (anti)periodic and /i = 
for free. The definition of B follows from (|^) 



B 







B 




and the definition of C follows from 



^ ^sinfc^(7^ ® 1) 



-C 




(24) 



(25) 



In this notation, following Neuberger's construc- 
tion leads to the free SDWF transfer matrix 



T = 



B-^/a^ 



-B~^C 



-CtS-i a^iC^B-^C - B) 



(26) 



Since we chose to set the staggered mass m 
in (||), then B is strictly anti-Hermitian, so T 
is as well (This is different from Wilson DWF and 
gives some idea why solving the zero mode prob- 
lem in (^) simplifies when solving for the field 
two sites away). In this case taking the as— i-O 
limit in order to identify the Hamiltonian does not 
make sense since the doublers are re-introduced 
(see section ||). On the other hand if we do not 
set the staggered mass w — > ^ then the as — > 
limit does not re-introduce the doublers and can 
be taken. However this term breaks the exact 
U(l)xU(l) symmetry. An alternative is to add 
an ^ — mo term as in section ^ This does not 
break the U(l)xU(l) symmetry, allows the as^O 
limit to be taken without reintroducing the dou- 
blers, but breaks the shift by one lattice spacing 
symmetry. 



7. Alternative actions 

The SDWF action considered here is not 
unique. Better actions may be constructed us- 
ing improved fields in the same spirit as with 
staggered fermions (see |19pG] and references 



therein). But even in the spin/flavor basis con- 
sidered here one could add the domain wall part 
in slightly diff'erent ways. For example one could 
have added to the standard staggered action the 
exact same term as the one added in Wilson DWF 
but multiplied by ^^^(ts ® Csa*)- 

8. Conclusions 

Staggered domain wall fermions (SDWF) have 
been constructed for the free theory. They de- 
scribe four flavors with exact U(l)xU(l) flavor 
chiral symmetry. The fuU SU(4)xSU(4) flavor 
chiral symmetry is recovered as the size of the 
extra dimension is increased. The addition of in- 
teractions for QCD is straight forward but it re- 
mains to be seen how well the SDWF properties 
are maintained at the stronger couplings where 
numerical simulations are done. It is hoped that 
SDWF may be useful for numerical simulations 
of QCD thermodynamics. 
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